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Abstract 

We estimated seedling, sapling, and mature tree densities for two historically defined grazed-classes 
(low and hi) of the Mad River Ranger District, Six Rivers National Forest, CA. Using Quercus garryana 
Hook, vegetation coverages on a geographic information system and a variable probability sampling 
scheme, 9 oak stands were randomly selected from each grazed-class. We employed a second sampling 
stage by selecting 3 simple random samples of 100-m 2 quadrants from each of the 18 oak stands selected 
in the first sampling stage. We found greater seedling and mature tree densities for the hi grazed-class 
and greater sapling densities for the low grazed-class. Sapling densities were roughly double and seedling 
densities about 5 times mature tree densities regardless of grazed-class. We suggest that increased grazing 
intensity creates favorable environments for seedling survival, but may ultimately reduce the number of 
seedlings transitioning to the sapling size-class. Our results showing roughly 2:1 sapling to mature tree 
ratios indicate that Q. garryana regeneration is occurring on these rangelands. 


The historical clearing of oak woodlands to max¬ 
imize forage yield for livestock grazing or to con¬ 
vert to agricultural production has been well doc¬ 
umented (Allen-Diaz and Holzman 1991; Bartolo- 
me et al. 1986; Bartolome and Standiford 1992; 
Muick and Bartolome 1987). Along with suburban 
encroachment, these land-use practices have been 
implicated as factors in the postulated decline of 
Quercus douglasii Hook. & Am (blue oak), Q. lob- 
ata Nee (valley oak), and Q. engelmannii E. Greene 
Engelmann oak regeneration (Griffin et al. 1987; 
Muick and Bartolome 1985; Pavlik et al. 1991). 
Regeneration problems manifested as bimodal age 
structure biased towards seedlings and mature trees 
with a lack of sapling age-classes have been cited 
for each of the above species (Bartolome et al. 
1986; Griffin 1971; Griffin 1976; Lathrop et al. 
1990). 

Quercus garryana Hook, has the longest north- 
south distribution among western Quercus spp. The 
northern range of Q. garryana extends onto Van¬ 
couver Island in Canada at 50°N latitude and the 
southern range runs into Los Angeles County at 
34°N latitude (Stein 1990). The coastal range of Q. 
garryana terminates in Marin County, but it is 
found further south along the western slope of the 
Sierra Nevada in a limited and disjunct distribution. 


The relationships between Northcoast Q. garry¬ 
ana populations and livestock grazing, land-use 
practices, and suburban encroachment have re¬ 
ceived limited attention (Griffin et al. 1987; Muick 
and Bartolome 1985; Reed and Sugihara 1986). Su- 
gihara et al. (1987) defined plant community types 
for Northcoast Q. garryana- dominated landscapes 
and Saenz and Sawyer (1986) studied grazing ef¬ 
fects on Q. garryana understory composition. Lit¬ 
erature treating potential domestic livestock grazing 
effects on Q. garryana is limited to a statewide sur¬ 
vey conducted by Muick and Bartolome (1987) and 
modeling efforts by Anderson and Pasquinelli 
(1984). Muick and Bartolome (1987) found no sig¬ 
nificant Q. garryana distribution patterns related to 
livestock grazing on the Northcoast. They assessed 
the adequacy of the existing Q. garryana sapling 
populations to replace trees lost through mortality 
by examining sapling to mature tree ratios on 13 
Northcoast plots. Sapling to adult tree ratios were 
less than 1:1 for Q. garryana in their study, yet 
75% of all individuals sampled were seedlings in¬ 
dicating no lack of seedling establishment, but a 
failure of seedlings to survive the transition to sap¬ 
ling size-classes. 

Bolsinger (1988) assessed the age structure of Q. 
garryana by drawing inferences from seedling and 
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Fig. 1 . Study sites located within Mad River Ranger Dis¬ 
trict, Six Rivers National Forest, Humboldt County, CA. 
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sapling counts from systematic plots established 
across California. He found 45% of Q. garryana- 
type plots contained no seedlings and 38% con¬ 
tained no saplings suggesting the potential for in¬ 
sufficient replacement populations or that present 
spatial distributions may be shifting. Anderson and 
Pasquinelli (1984) predicted declining Q. garryana 
distribution under “present trend” conditions. Their 
model parameterized livestock grazing, wildlife 
pressure, hre frequency, and weather conditions 
specific to the Northcoast’s Sonoma and Mendoci¬ 
no Counties. 

We designed this study to describe Q. garryana 
stand structure and to examine livestock grazing in¬ 
tensity effects on this structure. Specifically, we 
sought to determine whether differentially grazed 
stands contained significantly different densities of 
seedlings, saplings and/or mature trees. This infor¬ 
mation should provide insight about the regenera¬ 
tion status of Q. garryana on California’s North- 
coast. 

Methods 

Study site. We conducted this study June and July 
1995 on the Mad River Ranger District (MRRD), 
Six Rivers National Forest (SRNF), and Humboldt 
County, CA (Fig. 1). SRNF lies within sections of 
Humboldt, Del Norte, Siskiyou, and Trinity Coun¬ 
ties in the northern California coast range. MRRD 
is typified by a Mediterranean climate of cool, 
moist winters and warm, dry summers. Average 
January temperatures range from —2 to 4°C (28 to 
40°F) and average July temperatures range from 17 
to 27°C (64 to 80°F) (Oakeshott 1978). Annual pre¬ 
cipitation ranges from 127 to 152 cm (50 to 60 in). 


Table 1. Summary of Historical Grazing Allotment 
Data Used for Assigning Six Grazing Allotments to 
one of Two Grazed-Classes. 1 AU = Animal unit = 1 
individual. 


Allotment 

Mean 

annual 

1 AU 

Oak ha 

Oak 

ha-AU“ ] 

Grazed- 

class 

Long Ridge 

88 

1859 

21 

LOW 

Norris-Green 

108 

1839 

17 

LOW 

Barry Creek 

109 

1484 

14 

LOW 

Soldier Creek 

66 

616 

9 

HI 

Buck Mountain 

91 

766 

8 

HI 

Van Duzen 

166 

1166 

7 

HI 


Eighty percent of this precipitation occurs between 
November and April. Miles (1993) described 
MRRD soils as derived from the Franciscan Me¬ 
lange Complex with Quercus garryana stands oc¬ 
curring predominantly on the Oxalis-Hecker-Doty 
association. These are well-drained, pale brown 
loams with moderate to strong blocky structure and 
a slightly acidic nature. 

Grazed-classes. As an initial grazing allotment 
screening, we asked MRRD range personnel to 
identify 3 MRRD grazing allotments with high his¬ 
torical grazing intensity and 3 with relatively low 
historical grazing intensity. These allotments were 
to represent each of two grazed-classes (low or hi; 
Table 1). MRRD range personnel stated that each 
of the 3 low grazed-class allotments were histori¬ 
cally grazed very lightly, and had no permitted 
grazing activity for 5 to 10 years prior to this study. 
MRRD range personnel also confirmed that histor¬ 
ical permitted use (1950 to 1981) continued on the 
3 hi grazed-class allotments after 1981. 

We derived our grazed-classes by defining graz¬ 
ing intensities that were based on allotment-wise Q. 
garryana hectares (estimated using 1990’s GIS 
technology) because we believed that the majority 
of available forage, hence grazing pressure, oc¬ 
curred on the oak woodland/savanna vegetation 
type. We used annual animal units records (1950 to 
1981) and total Q. garryana hectare estimates for 
each allotment to quantify historical grazing inten¬ 
sities. Actual allotment-wise stocking rates were set 
by MRRD range managers using their assessments 
of suitable grazing hectares. Suitable hectares were 
comprised of oak savanna; coniferous forest; and 
small, open grassland areas known as “glades”. 
However, Sawyer et al. (1977) described the un¬ 
derstory layer of the Northcoast Douglas-fir—tan- 
oak type (Pseudotsuga menziesii (Mirbel) Franco— 
Lithocarpus densiflorus (Hook. & Am.) Rehder) as 
sparse, less developed, and dominated by shrubby 
L. densiflorus while others have confirmed this type 
of coniferous forest provides much less herbaceous 
productivity and available forage than oak wood¬ 
lands (Sharrow and Leininger 1982). While glades 
provide excellent forage, the fraction of the allot- 
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ments in this type is very small. Therefore, we as¬ 
sumed that oak woodlands provide the bulk of for¬ 
age for grazing animals and only used Q. garryana 
hectares in grazing intensity quantification. 

Sampling. Assuming that the total number of 
seedlings, saplings, and mature trees were a func¬ 
tion of the area of any oak stand, we employed a 
2-stage sampling design consisting of a variable 
probability selection scheme as the 1 st -stage and 
simple random sampling (SRS) as the 2 nd -stage. 
Nine Q. garryana stands were selected from each 
of 2 grazed-class subsets with probability propor¬ 
tional to the area of the remaining Q. garryana 
stands within that particular grazed-class. Hence, 
larger Q. garryana stands had higher selection 
probabilities. We incorporated these probabilities, 
known as l st -orderUnclusion probabilities (Overton 
and Stehman 1905), into Horvitz-Thompson (HT) 
estimation of the size-class densities for each 
grazed-class (Horvitz and Thompson 1952). We 
calculated sample-based sampling variance esti¬ 
mates for the two-stage HT estimator using the 
Sen-Yates-Grundy (SYG) method (Yates and Grun¬ 
dy 1953). This method requires that the 2 nd -order 
inclusion probabilities for the selected 1 st -stage 
units ( Q . garryana stands) be known. Second-order 
inclusion probabilities are the probabilities that any 
2 particular Q. garryana stands were both chosen 
during the selection process. We estimated 2 nd -order 
inclusion probabilities using the list sequential 
method outlined by Sunter (1977). The 2 nd sam¬ 
pling stage entailed SRS of three 100-m 2 quadrats 
from each unit (18 oak polygons) selected in the 1 st 
stage. We then scaled seedling, sapling, and mature 
tree estimates from each grazed-class, deriving 
mean values (per 100 m 2 ) for standardization and 
comparison ease. Ninety-five percent confidence in¬ 
tervals were constructed around each size-class 
density estimate for statistical comparisons. 

We obtained Universal Transverse Mercator 
(UTM) coordinates for each point selected during 
the 2 nd -stage of sampling from a geographic infor¬ 
mation system database. We arbitrarily decided that 
these UTM coordinates would serve as the south¬ 
eastern corner of each 100-m 2 sampling quadrat. 
Quadrats were located in the field with a handheld 
geographic positioning system. We delineated 
boundaries of each quadrat with a handheld com¬ 
pass and cloth tape. Seedlings (<1 cm basal di¬ 
ameter), saplings (<10 cm and >1 cm basal di¬ 
ameter), and mature trees (>10 cm basal diameter) 
were enumerated for each quadrat. Definition of 
size-classes followed Muick and Bartolome (1987). 
We did not age individual seedlings, saplings, or 
mature trees. 

We estimated several site variables at each quad¬ 
rat to determine potential environmental disparities 
among the two grazed-classes. These variables in¬ 
cluded 1) slope, 2) azimuth, 3) canopy cover, and, 
4) dominant herbaceous understory type. The sine 


of slope in degrees and the cosine of azimuth in 
degrees from N were multiplied to create a solar 
insolation indicator variable northness (sensu 
Borchert et al. 1989) which potentially varied from 
+ 100 to —100. Greater positive values of this vari¬ 
able indicate a N-facing steep slope while higher 
negative values indicate a S-facing steep slope. We 
ocularly estimated canopy cover at each site. We 
characterized the herbaceous understory type at 
each quadrat by making 6 randomly located her¬ 
baceous species cover estimates within the quadrat. 
Each estimate was made with a 10-point frame us¬ 
ing the first foliar intercept criterion for each pin 
placement (Heady and Rader 1958). We then clas¬ 
sified species into one of the following functional 
groups: 1) annual grass, 2) annual forb, 3) perennial 
grass, 4) perennial forb, or 5) woody perennial. 
Bare ground and dry organic matter 10-point frame 
hits were also recorded and categorized. 

We estimated edaphic characteristics at a ran¬ 
domly chosen 1 of 3 quadrats at each of the 18 Q. 
garryana stands. Soil pits were excavated for esti¬ 
mation of the following variables: 1) solum depth, 
2) rooting depth, 3) clay percent (by feel), 4) rock 
fragment percent (2-mm mesh sieve), 5) pH 
(LaMotte colorimetric method), 6) moist Munsell 
color value, 7) textural-class (by feel), and 8) avail¬ 
able water-holding capacity (USDA 1993). Para¬ 
metric tests were not performed to assess signifi¬ 
cant differences for these variables between grazed- 
classes, nor to correlate them to Q. garryana size- 
class estimates because they were collected with the 
variable probability sampling scheme described 
above. By definition, unequal probability samples 
do not meet the parametric test requirement of in¬ 
dependent, identical distributions between groups 
(Samdal et al. 1992). Therefore, we present contin¬ 
uous data in tabular form (Table 2) and discrete 
data in graphical (Fig. 2) form for grazed-class 
comparisons. Where tabular and graphical exami¬ 
nation suggested differences between the grazed- 
classes in a given variable, we tested joint bivariate 
distributions with a two-dimensional Kolmogorov - 
Smimov (2DKS) nonparametric procedure (Garvey 
et al. 1998). Two-DKS uses a bootstrapping tech¬ 
nique (500 Monte Carlo simulations) to iteratively 
compare randomly derived expectation matrices 
with an observed joint-distribution matrix. This 
technique is effective at assessing both linear and 
non-linear relationships among bivariate distribu¬ 
tions and makes no assumptions about functional 
responses. 

Results 

Comparison of 95% confidence intervals showed 
that seedling and mature tree densities were greater 
in the hi grazed-class while sapling densities were 
greater in the low grazed-class (Fig. 3). Mean es¬ 
timates are reported ± the 95% confidence interval 
and statistical significance determined by overlap or 
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Table 2. Sample Size (n). Means, and Standard Errors (SE) of Continuous Environmental Variables for Each 
Grazed-Class. 1 Units defined in text. 


Variable 

n 

HI 


LOW 

Mean 

SE 

Mean 

SE 

Canopy cover (%) 

27 

54.3 

4.5 

63.7 

4.4 

Slope (degrees) 

27 

17.5 

1.4 

26.9 

1.4 

Worthiness 

27 

2.6 

4.9 

-5.7 

4.6 

AWC 

9 

4.0 

0.5 

4.3 

0.3 

pH 

9 

6.2 

0.1 

6.2 

0.1 

Solum depth (cm) 

9 

115.7 

4.1 

118.4 

3.1 

Rooting depth (cm) 

9 

76.2 

2.0 

76.2 

0.7 

Clay (%) 

9 

15.4 

1.1 

17.4 

0.5 

Rock fragment (%) 

9 

46.0 

3.6 

42.7 

1.1 

Moist Munsell color value* 

9 

3.3 

0.1 

3.5 

0.1 


not. The seedling density estimate for the hi grazed- 
class was 33.5 ± 3.3 seedlings per 100 m 2 com¬ 
pared to 19.1 ± 1.2 seedlings per 100 m 2 for the 
low grazed-class. The mature size-class followed 
the same trend in that the hi mature tree estimate 
was greater than the low mature tree estimate (hi 
= 5.0 ± 0.5, low = 3.9 ± 0.4). Sapling size-class 
density estimates were the opposite. There were on 
average more saplings per 100 m 2 in the low 
grazed-class (10.8 ± 0.3) than in the hi grazed-class 
(9.3 ±1.1; Fig. 3). 

The greatest disparity in environmental variables 
between the two grazed-classes was average slope 
(hi = 17.5 ± 1.4, low = 26.9 ± 1.4; Table 2). 
Other continuous variables appeared quite similar 
so were not analyzed further than reporting means 
and standard errors. Examination of the dominant 
herbaceous understory type standard errors (Fig. 2) 
indicated similar distributions of this categorical 
variable between grazed-classes. Soil textural-class 
differences were evident; all 9 low grazed-class 
sites were classified as loams while 5 of 9 hi 
grazed-class soils were determined to be loams with 



matter 


1 clay loam and 3 sandy loams. Only 1 soil pit per 
site was dug, therefore no dispersion statistics were 
estimable from these data. 

Two-DKS results indicated that the joint distri¬ 
butions between slope and seedling density, slope 
and sapling density, and slope and mature tree den¬ 
sity (Fig. 4) were not significantly different from 
random (P = 0.60, 0.53, and 0.50; respectively). 

Discussion 

Greater seedling densities coupled with lower 
sapling densities in areas with greater grazing in¬ 
tensities suggest that herbivory of surrounding veg¬ 
etation may positively influence Q. garryana seed¬ 
ling survival and recruitment, but that incidental 
trampling and herbivory of the seedlings may dis¬ 
courage transition to the sapling stage. While this 
needs experimentation, there are experimental re¬ 
sults for Q. douglasii which corroborate our infer¬ 
ences and may help to explain our observations. 
Quercus douglasii and Q. garryana are taxonomi- 
cally and ecologically similar species. Both belong 


Grazed-class 

□ LOW 

□ HI 



Perennial Perennial Woody Bare 

grass forb perennial ground 


Functional group 


Fig. 2. Functional group cover means and standard errors (n = 9 for each grazed-class) as determined by 10-point 
frame estimates. 
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Fig. 3. Seedling, sapling, and mature tree density means and 95% confidence intervals for each grazed-class. 


to the white oak subgenus Lepidobalanus and oc¬ 
cupy apparently similar ecological sites—shallow, 
rocky hill slopes—in distinctive Californian cli¬ 
matic regions (Rundel 1979; Rundel 1986; Stein 
1990). Additionally, these species freely hybridize 
along climatic transition zones (Tucker 1979). 

In greenhouse experiments, Gordon et al. (1989) 
found greater Q. douglasii seedling emergence and 
growth responses in annual forb ( Erodium spp.) 
seeded plots compared to annual grass ( Bromus 
diandrus L.) seeded plots. They attributed this to 
increased rates of water stress in the annual grass 
plots. This illustrates that competition between Q. 
douglasii seedlings and neighboring herbaceous 
vegetation does occur, however, Gordon et al. 
(1989) obtained their results absent defoliation. 
Welker and Menke (1990) reported that defoliation 
of annual vegetation surrounding growing Q. doug¬ 
lasii seedlings was beneficial in reducing evapo- 
transpiration of competing vegetation, thereby re¬ 
ducing the rate at which the Q. douglasii seedlings 


encountered water stress. These results offer poten¬ 
tial explanations for our findings, whereby herbiv- 
ory is releasing Q. douglasii seedlings from re¬ 
source competition with surrounding vegetation. 
Resources being competed for are probably water, 
nutrients, and space. However, this phenomenon is 
likely overwhelmed at higher grazing intensities 
where grazing selectivity lessens and seedlings are 
depredated more or less incidentally. 

Hall et al. (1990) showed that at very high stock¬ 
ing rates, where animals are actually competing for 
space as well as forage, increased seedling depre¬ 
dation rates occurred. However, with reduced 
stocking rates. Hall et al. (1990) found little to no 
livestock preference for Q. douglasii seedlings. 
Given their similarities, these results likely hold for 
Q. garryana . Seedling depredation is probably a 
function of forage production and availability, ro¬ 
dent populations (Davis et al. 1990), and stocking 
rate. Experimental evidence for Q. garryana is 
needed. 



Slope (%) 


Slope (%) 


slope (%) 


Fig. 4. Joint bivariate distributions of slope and seedling density, sapling density, and mature tree density. 
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We did not measure microtopographic relief or 
estimate soil disturbance, but these factors may be 
important reasons for higher seedling densities. Be¬ 
cause livestock grazing creates soil disturbance 
(Kauffman and Krueger 1984), areas of higher 
grazing intensity tend to have greater micro-relief 
and therefore a higher probability of radicle inter¬ 
ception and penetration with the soil (Watt 1919). 

During autumn months, Q. garryana dormancy 
decreases interspecific competition for light and 
water between annual grasses and Q. garryana 
seedlings (Hibbs and Yoder 1993). This reduction 
in photosynthesis and growth allows seedling un¬ 
derground biomass persistence and promotes car¬ 
bohydrate storage in root systems (Hibbs and Yoder 
1993) while the annual grasses dominate available 
resources. Domestic stock grazing in the spring re¬ 
sults in annual grass defoliation as well as the in¬ 
cidental Q. garryana seedling defoliation and tram¬ 
pling. However, Hibbs, and Yoder (1993) have 
shown that Q. garryana is a prolific and adventi¬ 
tious sprouter. They found 20-year-old Q. garryana 
seedling roots with aboveground shoots less than 3 
years in age. Hence, should the Q . garryana seed¬ 
ling survive defoliation and water stress until an¬ 
nual grasses have been grazed or completed their 
life cycle, the Q. garryana seedling taproot and fine 
root system can eventually access resources at low¬ 
er soil depths where annual grasses are not com¬ 
peting. Indeed, we usually found a dense network 
of fine Q. garryana roots dispersed throughout the 
solum to a depth of 75+ cm (Table 2), while annual 
grass roots were not observed below the surface 10- 
cm. Seedlings able to establish these deeper roots 
should increase their survival chances via resource 
partitioning ( sensu Brown 1998). However, re¬ 
source partitioning is likely not achieved until 
Quercus seedlings have successfully negotiated in¬ 
terspecific resource competition mediated by cattle 
dietary preferences for herbaceous grasses and 
forbs. Coupled with an ability to survive many sea¬ 
sons of foliage removal via sprouting, resource par¬ 
titioning may account for greater Quercus seedling 
densities in areas with higher grazing intensities. 

Litter decomposition rate is another potentially 
important factor that may be differentially treated 
across grazing regimes. Desiccation is a major 
cause of germination failure in California’s pre¬ 
dominantly dry oak woodlands (Borchert et al. 
1989). However, the consistently moist, but cool 
Northcoast environs are more likely to result in low 
decomposition rates leading to higher litter accu¬ 
mulation with a concomitant increase in acorn rot. 
Opening of the herbaceous layer to increased solar 
insolation via phytomass removal by livestock may 
reduce the potential for acorn rot, thereby increas¬ 
ing germination and seedling success. These hy¬ 
potheses all need in situ experimental testing in Q. 
garryana stands. 

Significantly greater mature tree densities in the 
hi grazed-class may be the best explanation for both 


greater seedling densities and lesser sapling densi¬ 
ties in this group. It follows that more trees lead to 
more acorns leading to more seedlings. Transition 
from seedling to sapling is probably inhibited by 
lack of sufficient light and other resources due to 
crowding by mature trees combined with seedling 
defoliation and depredation by livestock. Allen- 
Diaz and Bartolome (1990) concluded that while 
recruitment of Q. douglasii seedlings was frequent, 
none of their seedlings made the transition into sap¬ 
ling size-classes. 

Why was mature tree density greater in the hi 
grazed-class? The Hi-grazed allotment group was 
generally found on gentler slopes that may have 
had some indirect effects on sapling to mature tree 
transition. For example, fires are known to bum less 
intensively on gentler slopes (Albini 1976) which 
might have provided a more favorable environment 
for tree development. Alternatively, deeper soils on 
sites with gentler slopes may have provided a mi¬ 
croclimate more conducive to sapling to tree tran¬ 
sition. However, our soil data revealed no evidence 
for this. Examination of maps and field observation 
revealed no obvious geomorphological or topo¬ 
graphical patterns that might account for these dif¬ 
ferences. 

Aside from the grazed-class comparisons of size- 
class densities, our results indicate that natural Q. 
garryana regeneration on these MRRD sites is not 
at risk. Representing roughly double the mature tree 
stock, densities of about 10 saplings per 100 m 2 
seems a surplus. Should one tree from a well- 
stocked stand die at any given time, several sap¬ 
lings should be available for replacement. To truly 
assess regeneration status, mortality and survival 
rates between each size-class must be known. Our 
single season study did not permit estimation of 
mortality rates even for seedlings. Our point-in¬ 
time estimates of age structure do provide compel¬ 
ling evidence that regeneration of this species is 
occurring. Annual assessment of seedling mortality 
rates coupled with creative techniques for estimat¬ 
ing sapling to tree transition probabilities and tree 
mortality rates are needed to verify our results and 
quantify the regeneration status of this resource. 

Conclusions 

Higher seedling but lower sapling densities with 
increased grazing intensity indicates that grazing at 
higher stocking rates does not affect Q. garryana 
seedling recruitment. There seems to be plenty of 
seedlings beneath and around tree canopies, but 
seedlings in the higher grazed areas are making the 
transition to the sapling size-class in lesser numbers 
than those in the more lightly grazed areas. How¬ 
ever, our results indicate that Q . garryana regen¬ 
eration is potentially occurring regardless of graz¬ 
ing pressures as evidenced by the roughly 2:1 sap¬ 
ling to tree ratio estimates found for both grazed 
classes. 
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